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We demonstrate the parametric amplification of the stochastic gravitational wave 
background during inflationary reheating and during axion/moduli oscillations. 
This enhances the detectability of the string/inflationary gravity wave signal, leav- 
ing a finger-print on the spectrum which might be found with future gravitational 
wave detectors. 



1 Introduction 

The nature of the cosmological background of gravitational waves is a subject 
of great interest at present, with many possible sources ranging from quantum 
fluctuations in a string-dominated phase Era or inflation U, to oscillations of 
cosmic string loops □. 

Here we focus instead on a powerful mechanism for distortion and am- 
plification of any existing gravitational wave background - namely damped 
parametric resonance due to oscillatory phases that the universe may -have 
undergone. Examples are provided by reheating at the end of inflation an 
oscillatory dilaton phase, or during coherent axion or moduli oscillations if 
they form a significant portion of the dark matter. 

We will discuss the amplification within the gauge-invariant Bardeen for- 
malism. Using the gauge-invariant and covariant electric and magnetic parts 
of the Weyl tensor gives similar results U. The evolution of the transverse- 
traceless (TT) metric perturbations h a b are naturally described by the Fourier 
mode functions h 6t f-, where e = {+, x} are the polarisation states. The hk f] 
satisfy: 

a ■ k 2 

h k + 3-h k + — h k = 0. (1) 

a a z 

Here a(t) is the scale factor of the universe which obeys the Friedmann Eq.: 

:)"-i"-5(5* + H- (2) 



"From now on we surpress the polarisation label. We also restrict ourselves to the case of 
flat spatial sections. 
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where k — 8nG and p is the relativistic energy density which we have specified 
to be in the form of a scalar field <fi, with potential V(4>). This gives us enough 
freedom to model both reheating physics and the oscillations of the axion 
condensate. Using Eq. (g), we can rewrite Eq. (|l]) asQ: 

d 2 (a 3 / 2 h k ) (k 2 , \ , , /2l , 
where p = 4> 2 /2 — V{4>) is the pressure. 



2 Parametric amplification of gravitational waves 

Now to illustrate parametric resonance, consider the n-dimensional first order 
system: 

y = P(t)y (4) 

where P is any matrix with period T. Then Eq. (|j) has n linearly independent 
normal solutions of the form: 



yi = Pi(t)e^ (5) 

where the fii are the characteristic/Floquet exponents of the system and the 
Pi are functions of period T. Then the n characteristic numbers defined by 
Pi = e<- iiT satisfy: 

pip2—p n = ex P {^J TrP{s)ds\ (6) 

with repeated characteristic numbers counted accordingly. The trace of P(t) 
is thus the crucial factor determining the existence of exponentially amplified 
modes. If TrP(t) > then eq. @ implies that: 

n?=ipi > i (7) 

which implies that at least one of the pi > 1 Pi > and hence by eq. (0) 

there is at least one unbounded, exponentially growing solution. 

The archetypal example is provided by the Mathieu equation (n = 2): 

y +[A~2q cos(2t)]y = (8) 

Let p = max{/ii, p%\. Then p is non-trivially related to the parameters (A, q) 
which span an instability chart consisting of an infinite hierarchy of resonance 
bands where p > OB. Defining e = AJ (2q) — 1, we plot p vs. (e, q) in Fig. (1) 
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using the piecewise quadratic approximation Em. This shows that for e < the 
resonance is particularly strong. 

The key point we wish to make is that Eq. (||) takes the form of the Math- 
ieu equation (||) when V((f>) = m^cf) 2 /2, as appropriate for studying chaotic 
inflationi. In this case, <fi evolves as ^>s'm(m^t). Then the pressure is: 

o 

777 

P= -^$ 2 cos(2m^) (9) 

yielding a Mathieu Eq. with parameters: 

fc 2 _ 3k$ 2 32k 2 

showing that in this case, unlike in the case of standard reheating with a posi- 
tive coupling constantly, e < is possible and gravitational wave amplification 
can be significant if $ ~ M v \. 

The effect of the expansion of the universe decreases $ and redshifts k, 
causing a decrease of both A and q, though e remains roughly constant. The 
decrease of q to below unity is particularly important in stopping the resonance, 
and there is thus a competition between the damping effect of the expansion, 
and the amplification due to resonance. 

One final important point regards the validity of temporal averaging. With 
an oscillating scalar field it is common to replace fi and p with their time- 
averages over an oscillation: ~p and p. In the case that V(<p) oc </> 2 , the average 
equation of state is that of dust, p = 0. From Eq. (||) this falsely predicts 
that there is no resonant amplification of the stochastic gravitational wave 
background during reheating or axion oscillations. This shows that temporal 
averaging is invalid. 

3 Applications 

3.1 Inflationary reheating 

The gravitational wave spectrum produced during inflation is nearly scale- 
invariant. However, during reheating via a second order phase transition this 
scale- invariance is broken and the rms value of the spectrum is amplified. From 
the form of q in Eq. (|Io|), this amplification is clearly strongly dependent on 
$, the initial amplitude of inflaton oscillations. This implies that the breaking 
of scale-invariance during reheating is much stronger in chaotic inflationary 
models than in new inflation, since $ is much larger in the former case. More 
discussion can be found inl9. 
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Figure 1: The Floquet index fj, on the stability- instability chart for Eq. (8). Notice the rapid 
increase of fi for decreasing e and increasing q. 



3.2 The axion and massive moduli 

The axio E is an oscillating scalar field and a natural cold dark matter can- 
didate. Unlike reheating, which lasts a very short time, the axion oscillations 
would last a large proportion of the universe's history, and hencejmight cause 
significant tensor amplification. The axion pontential is given byEoc3: 



V(4>) = A 4 



1 — COS 



fa 



(ii) 



with f a the axion decay constant and A = ,f a fn a , where m a is the axion mass. 
The standard QCD axion has A = A QCD ~ 200 MeV, f a ~ 10 12 GeV and 
gains a non-zero mass due to instanton effects at an energy around Aqcd^- 
There also exist massive moduli in supergravity and superstring theories with 
much less constrained parameters, for example one may take A ~ 10 16 GeV 
and f a ~ M p i Eil, with the moduli generically gaining mass at the epoch of 
supersymmetry breaking. 

To understand the implications of axion oscillations, let us approximate 
Eq. ( pi] ) by the first, quadratic, term in the Taylor series. We can then use 
the results of Eq. ( |l0| ) with the replacement — > A 4 // 2 so that roughly we 
have A ~ fc 2 / 2 /(A 4 a 2 ) and q oc $ 2 . For the values given above, this yields 



k 2 k 2 
hcD - 10 27 — , A moduU - 10" 26 ^ (12) 
a z a z 
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This implies that massive moduli are more likely to lead to large amplifications 
of the background gravity wave spectrum since e = A/(2q) — 1 < for a huge 
range of modes, while in the case of the QCD axion, only a tiny fraction of the 
modes, near k = 0, have negative e. 

On the other hand, only if the moduli or axions started with near-Planck 
expectation values, $ ~ M p i, will there be significant amplification in cither 
case. 



4 Conclusions 

We have shown that damped parametric resonance is important in understand- 
ing gravitational wave evolution during phases where a significant component 
of the energy density of the universe oscillates, such as during a second order 
phase transition or if the dark matter lies in an oscillating scalar field. 

This parametric resonance amplifies the resident stochastic background, 
changing the frequency dependence of the spectrum and enhancing the rms 
amplitude. This implies that the possibilities of detecting the stochastic back- 
ground of gravitational waves may be better than previously thought. In ad- 
dition there is the intruiging possibility of indirect detection of the axion or 
moduli via their finger-prints on the gravitational wave spectrum. 
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